Quantum oscillations near the metamagnetic transition in Sr3Ru207 
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We report detailed investigation of quantum oscillations in Sr3Ru207, observed inductively (the de Haas- 
van Alphen effect) and thermally (the magnetocaloric effect). Working at fields from 3 T to 18 T allowed us 
to straddle the metamagnetic transition region and probe the low- and high-field Fermi liquids. The observed 
frequencies are strongly field-dependent in the vicinity of the metamagnetic transition, and there is evidence for 
magnetic breakdown. We also present the results of a comprehensive rotation study. The most surprising result 
concerns the field dependence of the measured quasiparticle masses. Contrary to conclusions previously drawn 
by some of us as a result of a study performed with a much poorer signal to noise ratio, none of the five Fermi 
surface branches for which we have good field-dependent data gives evidence for a strong field dependence of 
the mass. The implications of these experimental findings are discussed. 



I. INTRODUCTION 

Quantum criticality is a topic of continued interest in 
condensed matter physics.' Quantum critical points (QCPs) 
have been found in various metallic systems, notably heavy 
fermion intermetallics and transition metal oxides, using tun- 
ing parameters such as hydrostatic pressure, chemical doping 
or magnetic field. They are usually accompanied by a signif- 
icant enhancement of the low temperature electronic specific 
heat, which is thought to stem from a quasiparticle mass renor- 
malisation due to strong quantum fluctuations associated with 
the continuous phase transition.^ 

QCPs have been studied by a large variety of methods, but 
surprisingly little is known about the microscopic nature of 
the processes involved. In particular, it is not clear if, in a 
multi-band system, the quasiparticles from all the bands will 
experience the mass renormalisation. The method of choice 
to determine electronic properties in k space has traditionally 
been the study of quantum oscillation measurements, which 
allow, in principle, the determination of the full Fermi surface 
and the associated quasiparticle massesj^ and can be studied 
over a wide range of magnetic fields. 

The subject of this paper, Sr3Ru207, gives an opportunity 
to probe the nature of quasiparticle mass enhancement in the 
vicinity of quantum criticality. It features a first order meta- 
magnetic transition which terminates in a critical point."* The 
depression of its critical temperature towards zero Kelvin us- 
ing a tilted magnetic field leads to a QCP, situated near 7.9 T 
when the field is aligned with the crystalline c-axis.^'^ Anoma- 
lous power laws in the resistivity were observed in its vicinity, 
typical of quantum critical systemsi^ A pronounced enhance- 
ment of the specific heat was reported near the QCP, with a 
logarithmic divergence of C/T at fields near the critical field. ^ 
Moreover, a previous de Haas-van Alphen (dHvA) study re- 
vealed the existence of at least five different electronic bands. 



of moderately high quasiparticle mass, two of which were re- 
ported to become enhanced when approaching the QCP with 
a magnetic field. ^ Away from the critical region, these masses 
were confirmed in an Angular Resolved Photoemission Spec- 
troscopy (ARPES) study by Tamai et al.^ 



In ultra-pure crystals, the approach to the QCP is signalled 
by a strong enhancement in the low temperature specific heat^ 
which is cut off by the formation of a new phase, within which 
transport properties develop a strong anisotropy under the ap- 
plication of tiny in-plane magnetic fields.' - In a recent let- 
ter we reported observation of dHvA inside this anomalous 
phase and demonstrated the importance of a chemical poten- 
tial shift in determining the field dependence of the dHvA 
frequencies »ii 

In this paper, we report a thorough investigation of quantum 
oscillations in Sr3Ru207 through the critical region. Perform- 
ing this experiment required the growth of a new generation 
of crystals, screened in a series of dHvA measurements in or- 
der to obtain the highest dHvA signal amplitude. Mean free 
paths of up to approximately 300 nm were measured, which 
enabled the detection of the de Haas-van Alphen effect down 
to fields as low as 2.4 T, as well as the observation of en- 
tropy oscillations using the magnetocaloric effect. We present 
here a detailed study of the quasiparticle masses of five out of 
six of the known bands as a function of the tuning parameter, 
the magnetic field. For each of these five bands, no signifi- 
cant increase in quasiparticle mass was detected. These results 
contrast with previous de Haas-van Alphen measurements by 
some of us, and so we explain the origin of this discrepancy in 
depth. We furthermore report measurements of the field and 
angle dependence of the measured dHvA frequencies. The 
implications of these new results are discussed. 
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FIG. 1: Schematic construction of the Fermi surface of Sr2Ru04, a 
and b, and of Sr3Ru207, c, d, e and /, where in the last two, the \/2 x 
\/2 reconstruction of the Brillouin zone (BZ) accompanied by back- 
folding was included, and the new BZ is represented by the dashed 
diamond. All axes were scaled by the parameters of the undistorted 
lattice. The red dots represent the location of a putative van Hove 
singularity. 



II. FERMI SURFACE 

The determination of the electronic structure in Sr3Ru207 
is not a simple problem, and has not yet been completely 
solved experimentally. It is a complex system partly due to 
the high degree of electronic correlations, and partly due to 
Fermi surface (FS) reconstructions arising from lattice distor- 
tions. To introduce the relevant issues, we first show the rela- 
tionship between the SrsRuiO? Fermi surface and the simpler 
one of its single layer relative Sr2Ru04. 

Sr2Ru04 possesses a body-centred tetragonal crystal struc- 
ture, of space group symmetry lA/mmm}'^ Consequently, it 
has fourfold rotational and inversion symmetry. Early work 
on Sr3Ru207 reported it to possess the same space group, 
but it was later found using neutron diffraction that this sym- 
metry is broken by a 7° rotation of the RuO octahedra.'^ '^ 
The resulting space group is Bbcb, corresponding to a B cen- 
tred orthorhombic structure with the a and b of sides being ef- 
fectively V2 larger than those of the undistorted crystal. This 



rotation of the octahedra has a significant effect on the elec- 
tronic system.—^ 

The FS in Sr3Ru207 can be derived schematically without 
involving band structure calculations, as shown in Fig.[T] One 
starts from the FS in Sr2Ru04, where three bands, a, p and 
y, are present. These originate from the hybridisation of the 
various d bands of the Ru atoms, the d^z, dy^ and d^y, which, 
due to crystal field splitting, are the ones crossing the Fermi 
level. The d^, and dy^ give rise to only weak dispersion ei- 
ther in the z and, respectively, y and x directions, due to quasi 
one-dimensional hopping. Consequently, the Fermi surface 
corresponding to these bands should be close to planar in the 
kxk, and kyk^ directions. These cross in certain regions of the 
Brillouin zone (BZ), where hybridisation gaps appear, and the 
sheets reconnect into closed surfaces and produce the a and 
P bands (Fig. [TJj.). The dxy orbital, however, allows hopping 
in both the x and y directions, and the corresponding Fermi 
surface is close to a perfect circle in the k^ky plane, the y band. 
The resulting FS with all bands is shown in Fig. This con- 
struction is consistent with the fourfold rotation symmetry of 
the lA/mmm space group. 

In Sr3Ru207, one expects each band to duplicate, due to 
bilayer splitting^**. Consequently, the result is slightly differ- 
ent than for Sr2Ru04. One starts from six surfaces (Fig.[T]:.), 
four that originate from the dxz, dy, orbitals, which reconnect 
at the points of crossing, and one should obtain a result sim- 
ilar to that shown in Fig. [T}/. This is what the FS would be 
without the \/2 x \/2 reconstruction due to the octahedral ro- 
tation. But since the BZ reconstructs into a square of half the 
area, back-folding of the bands occurs, shown in Fig. [T^. In 
this plot, many band crossings appear, and the way in which 
the surfaces reconnect is complex. According to dHvA and 
ARPES experiments,^'^ six orbits are expected, shown in Fig. 
[r|f . These take the form of square and cross shaped hole pock- 
ets in the centre, originating from the dxz and dyz orbitals, two 
lens shaped electron pockets at the M point, resulting from 
hybridisation between dxy and dxz, dyz, and a small pocket 
near the X point. These were labelled, respectively, ai, (X2, 
yi, p and y2 (see Fig. [Tlf). Finally, there is a small circular 
pocket centred on F, labelled 5, that has no direct analogy in 
Sr2Ru04. It has strong Ru c/^2_y2 character, and results from 
the depression in energy of a band that is unfilled in Sr2Ru04 
due to the lower symmetry of Sr3Ru207^i^ 



III. EXPERIMENT 

Sr3Ru207 crystals were produced using a floating zone 
image furnace, using the method published by Perry and 
Maeno.il 10 crystal rods were grown, from which more than 
50 individual samples were cut. Each was nearly cylindri- 
cal, 2 mm in diameter and 3 mm long, with the crystalline 
c-axis along the principal axis. These were the subject of a 
series of characterisation experiments, which enabled us to 
select the best four to use in these experiments, according to 
two criteria: disorder, due to vacancies and impurities, and 
phase purity, affected by the appearance of intergrowth re- 
gions with « 7^ 2 consecutive RuO planes (SrRu03, Sr2Ru04, 
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FIG. 2: Top panel Typical second harmonic dHvA oscillations mea- 
sured with B| |c-axis, at a temperature of 55 mK (see Mercure et al. 
for a discussion of the thermometry in this system)." Bottom panel 
Fourier transform of the data in the top panel between 5 and 6.5 T 
(shown as a red square). Five frequencies were observed, at 0.15, 
0.43, 0.90, 1.78 and 4.13 kT. 



Sr4Ru30io, etc.). 

Low disorder crystals were selected in a two-step process 
in which preliminary residual resistivity measurements were 
followed by short dHvA experiments performed in an adia- 
batic demagnetisation refrigerator. This featured the use of 
first harmonic detection, a DC field between 6 and 7 T and 
a temperature of 250 mK. The samples selected for the main 
experiments were those showing the largest amplitude of the 
1.78 kT frequency (see TableHll. 

In order to determine impurity phase volume fractions, 
we used the magnetic properties of the competing phases: 
superconductivity below 1.5 K (n = 1, Sr2Ru04),'^ ferro- 
magnetism below 160 K (n = oo, SrRu03),~~~ ferromag- 
netism below 105 K (« = 3, Sr4Ru30io).^ Since Sr3Ru207 
is paramagnetic^^ and produces only a small signal in suscep- 
tibility and magnetisation experiments, these properties were 
easily detected quantitatively for small amounts using zero 
field AC susceptibility (superconductivity) and magnetisation 
measurements using a Quantum Design MPMS SQUID mag- 
netometer (ferromagnetism). The samples used in the quan- 
tum oscillation experiments featured residual resistivities typ- 
ically of less than 0.4 /jilcm, and impurity phase volume frac- 
tions of a few % or less for Sr2Ru04 and 0.1% for SrRu03 
and Sr4Ru30io. 

The de Haas-van Alphen effect was measured using a dilu- 
tion refrigerator located in the Cavendish Laboratory in Cam- 
bridge, equipped with an 18 T superconducting magnet. The 
samples were placed in well compensated coil pairs, each 
comprising 1000 turns of 25/jm diameter Cu wire. These were 
placed onto a rotation mechanism which allowed the measure- 



ment of the magnetic susceptibility at angles between the field 
and the crystalline c-axis of up to about 60°, both along the 
[100] and [110] directions. The coil pairs were connected to 
the input coils of low temperature transformers, which were 
themselves coupled with room temperature preamplifiers. The 
magnetic susceptibility was measured using first and second 
harmonic detection, with digital lock-in amplifiers and an ex- 
citation field of 20 G at 7.9 Hz. These values for the modu- 
lation field were chosen in order to ensure that eddy current 
heating in the samples was negligible. They were selected 
following an extensive study which demonstrated both that 
the measured masses were independent of modulation fields 
up to this level, and that, under the conditions selected for 
the experiments on Sr3Ru207, the well-known masses of all 
sheets of Sr2Ru04 were obtained correctly. The noise floor at 
the compensated coils was of around 30 pVrms/\/Hz, when 
using an analog magnet discharge field sweep method, and 
150 pVrms/\/Hz when using linear field sweeps. 

For the majority of the measurements presented here, sec- 
ond harmonic detection was employed in order to minimise 
the effects of the large magnetic background near the meta- 
magnetic transition. A series of sweeps was performed be- 
tween 18 T and 2.5 T, at rates of approximately 0.04 T/min 
for B> 10 T and 0.02 T/min for lower fields. These were car- 
ried out at 10 different temperatures between 90 and 500 mK 
with the field aligned with the crystalline c-axis, and then at 
angles from -15 to 55° using a step of about 1.6°, at a temper- 
ature of approximately 55 mK. For one aspect of the data set, 
the field dependence of very low dHvA frequencies at low 
fields, we performed additional measurements. Oscillations 
were studied in St Andrews using the magnetocaloric effect in 
which the sample temperature is monitored while the field is 
swept. The relaxation times associated with this kind of mea- 
surement favour low frequencies over higher ones. To follow 
those up with dHvA we worked with first harmonic detection, 
maximising the low frequency signal at the expense of a large 
magnetic background which restricted the experiment to fields 
below 7 T. 



IV. DHVA DATA, SPECTRA AND QUASIPARTICLE 
MASSES 

In Fig. |2]we present susceptibility data between 3 T and 
18 T and the dHvA spectrum of data between 5 T and 6.5 T. 
The top panel shows that dHvA oscillations can be detected 
throughout the available field range. Three metamagnetic sig- 
natures can be seen near 8 T, which consist of asymmetric 
peaks centred at 7.5, 7.9 and 8.1 T. The first corresponds to 
a metamagnetic cross-over and the second to two first or- 
der transitions, as determined using the out-of-phase com- 
ponent of the first harmonic AC susceptibility (not shown). - 
The region between 7.9 and 8.1 T coiTesponds to the anoma- 
lous phase discussed by Borzi et alJ^ and our recent work 
in which we reported dHvA oscillations."'^^ The dHvA spec- 
trum shown in the lower panel of Fig. |2] features five peaks, 
situated at 0.15, 0.43, 0.90, 1.78 and 4.13 kT. In previous 
ARPES work by some of us,^ we associated these with bands 
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TABLE L Electronic data for all six bands observed in Sr3Ru207. 
The type refers to whether bands are electron-like (e) or hole-like 
Qi) and the multiplicity refers to possible bilayer splitting times its 
number of occurrences inside the BZ. F is the dHvA frequency, the 
number of electrons refers to the fraction of the BZ area covered by 
a band, which includes the multiplicity. Finally, m* jnie is the quasi- 
particle mass and Yf/, is the associated contribution to the electronic 
specific heat coefficient, including the multiplicity. 



labelled respectively p, 5, yi, ai and a2, where the first three 
correspond to electron orbits and the last two to hole orbits. 
A sixth frequency also exists at 110 T, as will be discussed in 
sec. |Vl Note that the splitting of the peak at 0.9 kT into four 
peaks is associated with bilayer splitting of this band, as well 
as to a normal corrugation of this branch of the FS, while the 
splitting of the peak at 4. 13 kT is due to magnetic breakdown, 
as detailed in section VI. 

Taking into consideration all the bands reported, the 
ARPES map of the BZ and the predicted multiplicity and 
bilayer splitting of these bands, we were able to put for- 
ward a model which allowed us to perform the Luttinger 
sum and account for the total electronic specific heat coeffi- 
cient in Sr3Ru207, which was reported by Ikeda et al. as 1 10 
mJ/molRuK^."^ It involves assigning multiplicities to bands 
and determining whether they are electron or hole-like. Table 
J] shows a summary of this model, along with the dHvA fre- 
quencies of the six bands, the associated number of electrons 
which occupy these, their associated quasiparticle masses and 
their contribution to the total specific heat. The number of 
electrons was calculated by normalising the dHvA frequencies 
by that corresponding to the in-plane area of the BZ, 13.7 kT, 
multiplied by the multiplicity of the bands, and taking into ac- 
count the spins. For hole-like bands, these are expressed as the 
number of filled zones minus the fraction of electrons covered 
by these bands. The contributions to the specific heat were 
calculated using the quasiparticle masses, times the multiplic- 
ity, and a conversion factor of 0.74 mJ/molRuK^i^ Moreover, 
a number of electrons are included in zones that become filled 
in the process of the back-folding of the bands, which ac- 
counts for 8 electrons. The total number of electrons counted 
in this way nearly matches the required amount, 15.71 com- 
pared to 16 electrons, while the calculated specific heat is con- 
sistent with that measured previously, 103 ± 8 compared to 
llOmJ/molRuK^. 
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FIG. 3: (a) The change of sample temperature during a sweep of 
the magnetic field at 0.04 T/min (blue) and -0.04 T/min (red) as a 
function of magnetic field /joH. The average sample temperature is 
160 mK. (h) Fourier transform in II^qH of the blue trace in panel 
(a). In green is given the Fourier transform of a dHvA measurement 
over a comparable magnetic field region at base temperature. The 
label (A) indicates the new frequency found at 110 T. (c) Tempera- 
ture dependence of the amplitude of the 3 kT frequency of Sr2Ru04 
as measured by magnetothermal oscillations. The red curve is a fit 
of the temperature derivative of the Lifshitz-Kosevich formula to the 
data, (d) Amplitude of the 110 T peak in the Fourier transform of 
the data in panel (a) as a function of sample temperature Ts- The red 
curve is a fit of the temperature derivative of the Lifshitz-Kosevich 
function to the data. This fit gives a mass OTnor for the newly ob- 
served frequency of 8(1) nie- 



V. IMAGNETOTHERMAL DETECTION OF LOW FIELD, 
HIGH MASS OSCILLATIONS 

One of the predictions of the ARPES work of Tamai et ali^ 
was the possible existence of tiny Fermi surface sheets associ- 
ated with extremely flat bands that they labelled 72. Observing 
tiny frequencies against the background of many other oscil- 
latory components is not easy, especially since the flatness of 
the bands means that for any given field range, the signal size 
is likely to be suppressed relative to that of other components 
of the overall signal. This effect is exacerbated when work- 
ing with second harmonic dHvA, in which relative signals for 
low frequencies go as F^. In order to search for the oscil- 
lations, we profited from bespoke equipment developed in St 
Andrews for the study of the magnetocaloric effect.^ In con- 
trast to second harmonic dHvA, the relative signal size will be 
independent of F in principle, but subject to a damping fac- 
tor that will grow rapidly with F since more rapid oscillations 
will be damped more strongly in a measurement system with 
long intrinsic relaxation times. 

In Fig. [3^ we show low frequency oscillations detected by 
the magnetocaloric effect for up- and down-sweeps of the field 
in red and blue respectively. The sign reversal of the tempera- 
ture change is a key diagnostic that these are genuine magne- 
tothermal oscillations. If, for example, they were due to sam- 
ple heating due to eddy currents subject to the Shubnikov-de 
Haas effect, they would have a sign independent of field sweep 
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FIG. 4: Field dependence of the dHvA spectra observed with the field 
aligned with the direction of the c-axis. The colour scale corresponds 
to the amplitude of the spectrum. 



direction. Their Fourier transformation is shown in the black 
trace in Fig. [S}?. The previously undiscovered peak at 110 T is 
shaded. In the same panel, the green trace shows the Fourier 
spectrum of a subsequent dHvA run taking data on the first 
harmonic with a modulation field of 3 gauss. While confirm- 
ing the existence of the HOT peak, it emphasises the relative 
signal enhancement that magnetocaloric oscillations can pro- 
vide for very low frequency signals. 

An important aspect of magnetocaloric oscillations is that 
the amplitude of the signal is proportional to dM/dT rather 
than to M directly so the temperature dependence expected 
of the oscillations is that of the first thermal derivative of the 
usual Lifshitz-Kosevich form (see equation [T] below). This is 
shown in red in Fig. [3}; which shows the result of a 'con- 
trol' experiment detecting magnetocaloric oscillations from 
the well-known a sheet of Sr2Ru04. In spite of the different 
temperature dependence of the raw signal, the mass deduced 
agrees well with that obtained by standard dHvA, providing 
a useful verification of the calibration of the magnetocaloric 
apparatus. 

In Fig. [3}l we show the temperature dependence of the mag- 
netocaloric signal of the new frequency at 110 T, along with 
the fit that allows an effective mass estimate of 8±1 m,,. A 
subsequent first harmonic dHvA study involving careful sig- 
nal averaging increased the precision of the measurement to 
7.2±0.2otj.. This mass coiTesponds to a Fermi velocity of only 
9 kms^', a factor of six smaller than that of the y sheet of 
Sr2Ru04. 



not constant along the field range, but change as the system 
crosses the metamagnetic transition. Variations are observed 
as one increases the magnetic field and are especially signifi- 
cant above 7 T. 

In Fig. |4]we present an analysis of the field dependence of 
the dHvA frequencies along the field range. It was calculated 
using multiple Fourier transforms taken over short field ranges 
of equal inverse field width. 100 such windows were used, of 
width of 0.02 T^', which overlap and are distributed evenly 
in inverse field space^'^. The intensity of each spectrum cal- 
culated is plotted on a colour scale. The hidden central area 
coiTesponds to the metamagnetic region, where a large non- 
oscillatory signal arises (see Fig.|2]near 8 T). We have recently 
succeeded in detecting a tiny oscillatory signal on top of the 
large background, as reported in ref. [11]. 

On the low field side of the metamagnetic transition, the 
same five dHvA frequencies are observed as in Fig. |2] (ex- 
cluding 72). However, near 7 T, increases are observed, par- 
ticularly with the 1.78 and 4.13 kT frequencies, which cor- 
respond to the ai and a2 bands. This is discussed in detail 
in ref. [11], and corresponds to a change in the size of these 
parts of the FS, due to a change in Fermi level. Note that the 
fact that the ai and a2 bands remain of similar size in high 
fields compared to low fields exclude any further V2 x \/2 
reconstruction of the FS at the metamagnetic transtion. 

One can also observe in Fig. |4] that the data become more 
complex as the magnetic field increases, and that most dHvA 
peaks split into several components. The new frequencies 
might originate from magnetic breakdown, activated as the 
field increases. In particular, for the case of a2 (see Fig.[l]f.), 
symmetry considerations allow four magnetic breakdown re- 
gions with the Yi pocket, which increase the cross-sectional 
area by small amounts. This leads to a large number of pos- 
sibilities, each new orbit occurring with a slightly different 
probability. In such a case, the amplitude of the dHvA signal 
becomes divided between a large number of possible frequen- 
cies all situated near 4 kT. The same reasoning holds with the 
tti pocket, from which the electrons can tunnel towards the 
a2 pocket and back in four different locations, leading to a 
similar number of possibilities. 

It is difficult to explore this in any more detail, since the 
topology of the Fermi surface obviously changes slightly 
across the metamagnetic transition and is not known above 
8 T. It is thus impossible to know the magnetic breakdown 
field values precisely. We note however that, as can be seen 
in Fig. m the frequency splitting that we associate with mag- 
netic breakdown appears at fields lower than that of the meta- 
magnetic transition, with the 4. 1 kT peak splitting into three 
or more between 6 T and the metamagnetic transition. This 
suggests that the breakdown and the metamagnetism are not 
directly related^. 



VI. FIELD DEPENDENCE OF DHVA FREQUENCIES AND 
POSSIBLE MAGNETIC BREAKDOWN 

As was already suggested by Borzi et al.^ and described 
in more detail in our recent work^!^ the dHvA frequencies are 



VII. ANGLE DEPENDENCE OF DHVA FREQUENCIES 

Measurements of the angle dependence of the dHvA sig- 
nal allow one, in principle, to determine the three dimensional 
topology of the FS. In quasi-two dimensional materials, or- 
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FIG. 5: Top panel Angular dependence of the dHvA spectra ob- 
served in the low field side of the metamagnetic transition, where 
the colour represents the spectral intensity, while the horizontal and 
vertical axes represent the angle and the dHvA frequency multiplied 
by cos(6) respectively. Fourier transforms were taken between 5 and 
6.5 T. Bottom panel Angular dependence of the spectra in the high 
field side of the metamagnetic transition. Fourier transforms were 
taken between 10 and 18 T. 



bits are open in the direction, and the frequencies follow 
a simple angular dependence of the form F{Q) = Fq/cosG, 
where Fq is the frequency when the field is aligned with the c- 
axis direction. This is observed in the case of Sr3Ru207, for 
fields lower than that of the metamagnetic transition. Fig. |5] 
top panel, presents Fourier transforms of the data between 5 
and 6.5 T measured at angles between and 54°, plotted with 
the frequency multiplied by cos(9), such that cylindrical FS 
pockets appear as horizontal lines. The amplitude of the spec- 
tra are represented in colour in vertical strips for each angle. 
All of the FS pockets observed follow almost perfectly the ex- 
pected behaviour for a 2D system, and small deviations are 
due to their small warping^'. This agrees well with the pre- 



dictions made by band structure calculations.'^' ° 

The situation is not quite the same at fields above that of 
the metamagnetic transition. We found that frequencies are 
not continuous functions of angle; some of them appear while 
others vanish at specific angles. In particular, the frequency at 
1.6 kT, related with the ai pocket, disappears between 5 and 
10°, while that at 1.0 kT, probably related to the yi pocket, is 
not present for B ||c-axis and appears at a similar angle. We 
observed a similar behaviour with the frequency near 0.43 kT, 
which appears near 20°. Finally, frequencies near 3 kT appear 
briefly at low angles. These phenomena may also be related 
to the magnetic breakdown effects mentioned earlier. Effec- 
tively, when a two dimensional FS possesses small corruga- 
tions, magnetic breakdown fields can vary as the field angle 
is changed, due to a changing ^-space distance between FS 
pockets. If such a distance decreases significantly, magnetic 
breakdown becomes overwhelmingly probable, to the extent 
that the original orbit becomes less probable, and may escape 
detection. Varying magnetic breakdown gaps might therefore 
make the angle dependence of dHvA frequencies appear as if 
the FS changed topology with angle, while no such changes 
occur in reality. 



VIII. QUASIPARTICLE MASSES 

The enhancement of the specific heat previously reported 
near the metamagnetic transition is expected, in a metal, to be 
reflected in the quasiparticle masses measured from quantum 
oscillations.^ Borzi et al. reported strong enhancements of the 
quasiparticle mass of the ai and a.2 bands, below and above 
the metamagnetic transition.^ It was moreover expected that, 
although not possible to determine at the time, the masses of 
the other bands should be enhanced as well, stemming from 
quantum critical fluctuations coupling to the whole Fermi sea. 
Indeed, the desire to confirm and extend the findings of ref. (01 
was the main factor motivating this new study using samples 
and detection methods designed to deliver an improved signal 
to noise ratio. To our surprise, we did not reproduce our pre- 
vious findings. In contrast, we demonstrate here that no such 
enhancement occurs for the five bands which we detected in 
second harmonic dHvA; ai, a2, 8, yi and p (see Fig.[Tl/'). We 
furthermore show that the data analysis method used by us in 
ref. 101 to calculate quasiparticle masses can lead to erroneous 
results if used with data of low signal to noise ratio, as we 
discuss in detail in appendixlAl 

In Fig. |6]we present the field dependence of the quasipar- 
ticle masses of the dHvA frequencies at, in order from top 
to bottom, 0.15, 0.43, 0.90, 1.78 and 4.13 kT, which corre- 
spond, respectively, to p, 5, yi , a.\ and Oa (excluding y2), us- 
ing a magnetic field aligned with the crystalline c-axis. The 
field dependent masses were calculated using 100 field win- 
dows of equal inverse field width of 0.01 T^^ between 4 and 
18 T, over which Fourier transforms were performed for all 
measured temperatures. The amplitude of the oscillations at 
specific frequencies was obtained by integrating peaks in the 
power spectra. These were subsequently analysed using the 
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FIG. 6: Quasiparticle mass of the five frequencies identified in the 
low field side, as function of applied magnetic field. These corre- 
spond respectively to the bands ai, ai, 5, Yi and (3. The anomalous 
phase region, situated between 7.9 and 8. 1 T, is indicated with verti- 
cal red lines. The fits were calculated over field ranges of 0.01 T^' 
except for the 0.15 kT frequency, for which a width of 0.02 T^' was 
necessary in order to properly resolve the corresponding peaks in the 
spectra. Gaps in the data are due to the amplitudes becoming too 
small to fit because of beat patterns in the data. 



standard Lifshitz-Kosevich (LK) relation 



LK{Bo,T)^A 



Cm*T/BQ 
sinh(COT*r/Bo)' 



Bo 



I/B1 + I/B2' 



(1) 



where A is a temperature independent factor, Bq is the inverse 
of the average inverse magnetic field, C contains all universal 
constants and is approximately equal to 14.7 T/Kmg and m* is 
the quasiparticle mass. Note that the change to Bo due to the 
increase in internal field around the metamagnetic transition is 
of less than 0.2% and was neglected. Non-linear fits of this re- 
lation to the data were performed using two parameters, A and 
m*, for each field window, and each fit was carefully inspected 
by eye, rejecting field windows which produced inappropriate 
data distributions. The resulting masses are plotted in Fig. |6] 
as a function of Bo- The error bars were obtained from the 
standard non-linear regression procedure. Fig. |7] shows four 
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FIG. 7: LK fits performed in various field ranges for three of the six 
detected frequencies, 0.15, 1.78 and 4. 13 kT, which correspond to the 
P, Ml and a2 bands. As in the data of Fig. |6l for the case of 0. 1 5 kT, 
inverse field windows of a width of 0.02 T^' were used, while for the 
other two the width was of 0.01 T-i. The field regions are indicated 
in colour boxes in the respective insets, and the same colour schemes 
are respected for the data. Integrated peaks in the Fourier spectra are 
represented with various symbols, while fitted LK curves are shown 
with solid lines. The values in the legends correspond to the extracted 
quasiparticle masses. 



of these LK fits for the frequencies at 0.15, 1.78 and 4.13 kT 
in selected field regions, which are shown as colour boxes in 
the respective insets. The fits are mostly very good. 

It is clear from the data in Fig. |6] and from the fits in Fig.|7] 
that no strong enhancements of the quasiparticle masses of the 
bands p, 8, yi, ai and a2 are present near the metamagnetic 
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transition. Changes of at most 4 electron masses were seen, 
which originate from fitting errors when the signal to noise 
decreases at, for instance, nodes of beating patterns. For the 
case of Y2, situated at 0. 1 1 kT, its resolution from the peak at 
0. 15 kT required a field range on the low field side of around 
4 to 7 T, making the calculation of the field dependence of 
its mass impossible. The second harmonic data shown here 
suppress its signal so much in relation to that of the 0.15 kT 
frequency that the temperature dependence of the latter could 
be successfully isolated. Note furthermore that the absence of 
data for the peaks at 0.43 and 0.90 kT on the high field side of 
the metamagnetic transition is due to the absence of signal in 
that region for these bands (see Fig. H]). The signal for these 
bands appears at higher angles (see Fig. |5]l, suggesting that 
another process suppresses the amplitude of their quantum os- 
cillations at low angle but that these bands are still present in 
that region. This result has been reproduced several times, and 
additional verifications have been performed using first har- 
monic detection and a lower modulation field (3 G compared 
to 20 G). 

The absence of a quasiparticle mass enhancement as a func- 
tion of magnetic field in Sr3Ru207 is a surprise to us. Both the 
electronic specific heat^ and the weight of the Fermi liquid 
term"* show approximately consistent enhancements that, to- 
gether, imply an average rise of the quasiparticle mass of a 
factor of between 2 and 3 between low fields and 7.5 T. It was 
therefore vital that we checked the sensitivity of our analy- 
sis techniques to an enhancement of this scale for all the de- 
tectable frequencies in our signal. We decided that the only 
way to test this in practice was using simulated dHvA data. 
Starting with the same parameters as those of our experiment 
(including modulation field, data sampling rate, relative am- 
plitude factors for the different oscillatory components and 
estimated sample mean free path), we modelled a composite 
signal comprising the P, 5, yi, ai and a.2 frequencies to which 
we added random uncorrelated Gaussian noise of a similar 
magnitude to that present in our experiments. 

The simulated data were subsequently analysed using ex- 
actly the same procedure as for the experimental data. The 
first data set used a field varying quasiparticle mass for the five 
bands simulated, using the experimental value away from the 
metamagnetic transition, and a strong, gradual enhancement 
shown by the red lines in Fig. [8] The onset of this increase was 
determined approximately from that of specific heat data.^ In 
contrast, the second simulation featured no quasiparticle mass 
enhancements, but instead an increase in scattering rate of a 
factor of 2 at 8 T in order to replicate the observed decrease in 
signal at the metamagnetic transition (see Fig. |2] top panel). 

In Fig.|8]we present the results of the analysis of the simu- 
lated data featuring mass enhancements for all the bands. The 
red curves indicate the values of the mass used in the calcu- 
lations. We observed that depending on the signal to noise 
of the various oscillatory components, the extracted masses 
track the original ones very well for mass changes of up to 
around 10 m^. for the component with lowest signal to noise 
ratio (0. 15 kT), to above 25 nie for that with the highest quality 
signal (1.78 kT), all of which are much larger than the largest 
mass changes measured in our data (of at most 4 nie). We thus 
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FIG. 8: Quasiparticle mass of five simulated dHvA frequencies sim- 
ilar to those detected in our experiments, calculated using the same 
method as that used in Fig.|6] The data was modelled using quasipar- 
ticle masses which were strongly enhanced near 8 T. The oscillations 
were calculated using the masses and mean free path extracted from 
the experiment, as well as amplitudes, field dependent frequencies 
and beat patterns similar to those in the data and uncorrelated gaus- 
sian distributed noise was added to the oscillations. 



conclude that if such mass changes had been present in our 
dHvA data, we would have detected them. 

In Fig.|9]we show the mass analysis of the second data set 
which does not feature mass enhancements, and yields the ap- 
propriate mass values. Appendix lAl describes how the previ- 
ous incorrect result probably arose as a consequence of using 
a free offset as a third parameter in the Lifshitz-Kosevich fit to 
attempt to process data with a low signal to noise ratio. Since 
the current data had a much better signal to noise ratio, the 
use of either two or three parameter fits did not produce sig- 
nificantly different results. 

In a third set of tests we experimented with assigning a di- 
verging mass individually to different frequencies, in case the 
entire specific heat change were due to a diverging mass on 
only one Fermi surface sheet. In doing so, we confirmed that 
if this were happening, our data analysis procedures would be 
sensitive to it. The only exception to this came when trying 
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FIG. 10: Quasiparticle mass extracted from the merged peak com- 
posed of the frequencies at 0. 1 1 and 0. 15 kT, originating respectively 
from the Y2 and (3 bands, where in the simulation, the mass of Y2 band 
is strongly enhanced while that of the P band is not. The mass en- 
hancement is washed out by the P band, which possesses a mass of 
5.5 nip. 



this Fermi surface sheet on the basis of the current data. 
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IX. CONCLUSION 



FIG. 9: Quasiparticle mass of five simulated dHvA frequencies sim- 
ilar to those detected in our experiments, calculated using the same 
method as that used in Fig.[6] The data was modelled using a constant 
quasiparticle mass and a scattering rate which was strongly enhanced 
near 7.9 T, driving the amplitude to zero at this field value, similar to 
the experimental observations. 



to include a diverging mass for a frequency at 0. 11 kT. Here, 
it is impossible to perform the field-dependent mass analysis 
and still resolve that component of the signal from the one at 
0.15 kT. In first harmonic dHvA, the integrated weight of the 
0. 1 1 kT component is much less than that of the 0. 15 kT com- 
ponent, and the effect of the 0.11 kT component experiencing 
a rapid change of amplitude with temperature is masked by the 
much weaker temperature dependence of the 0. 15 kT compo- 
nent. This is demonstrated in Fig. [TOl which shows the re- 
sult of a simulation incorporating the experimental oscillatory 
amplitudes in which the mass of 72 is enhanced by an order 
of magnitude between low / high fields and the transition re- 
gion. The relative integrated weights of our magnetocaloric 
data are similar at our lowest temperature of measurement, 
but for that apparatus the base temperature is only 100 mK, 
and resolving a strongly field dependent mass for the 0.11 kT 
component would require lower temperature data. We con- 
clude, therefore, that we can neither definitely detect nor rule 
out a strong field dependence of the quasiparticle mass from 



The first two pieces of new information that we have pre- 
sented in this paper concern the magnetic field and field angle 
dependence of the dHvA frequencies both below and above 
the metamagnetic transition fields in Sr3Ru207. Although the 
data contain some unusual features, it seems possible to ac- 
count for these in terms of previously known effects. Field- 
dependent frequencies are a consequence of back-projection 
in the presence of itinerant metamagnetism.^^ .'' The 'frag- 
mentation' of the observed frequencies at high fields can be 
qualitatively accounted for by magnetic breakdown effects, al- 
though we do not have sufficient information to perform the 
full quantitative analysis that would be required to prove that 
this is the sole cause of the observations. 

The finding that we observe essentially no systematic field 
dependence of the quasiparticle masses is more surprising. 
Taking a conventional viewpoint, it is tempting to assume that 
since the specific heat shows a large change over the same 
field range, some quasiparticle mass enhancement must be 
taking place, and hence that it must be associated with the 
72 frequency for which we have no reliable dHvA information 
on the field dependence of the mass. Adopting that assump- 
tion would imply an order of magnitude change of its mass 
between low field and 7.5 T^^. It would be dangerous to draw 
a firm conclusion on this point in the absence of firm data; this 
aspect of our work motivates further experimental and theoret- 
ical study. 
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Appendix A: Systematic errors in dHvA mass analysis 

We present in this appendix a discussion of dHvA quasi- 
particle mass analysis methods, and how certain assumptions 
in the analysis method can lead to large systematic errors. In 
particular, we try to explain how the analysis performed in 
the work of Borzi et alJ- led to the apparent observation of a 
strong quasiparticle mass enhancement of the ai and a2 bands 
(at 1.78 and 4.13 kT). 

The usual analysis method of the temperature dependence 
of the amplitude of dHvA oscillations consists of performing 
Fourier transforms of data for all temperatures over the same 
field range. Integrals are performed over a specific peak in 
the resulting power spectra, of which the square root is taken. 
The field dependent amplitudes are then analysed using a non- 
linear fit of the standard LK expression, of the form given 
in eq. [T] This function saturates at the lowest temperatures 
towards the parameter A of Eq. [T] and exponentially decreases 
at high temperatures. 

However, in quantum oscillation data which feature ap- 
preciable levels of noise compared to the signal, the inte- 
gral of the power spectra never decreases to zero amplitude 
at high temperatures. Effectively, if we write the signal as 
Acos{Fo/B + (^) +f{B), where f(B) consists of random noise, 
the power spectrum W{F) as a function of frequency F of the 
Fourier transformed signal can then be written as 



W{F) = \d{F-Fo)+f{F)\ 



(Al) 



where 5(F — Fq) corresponds to a peak centred at Fq and f{F) 
is the Fourier transform of f{B). At high temperatures, the 
peak vanishes and one is left with W{F) — \ f{F)\^, which is 
not zero. In the analysis of ref. [f] we subtracted the apparent 
resulting offset in the high temperature range of the tempera- 
ture dependence of the dHvA amplitude by using a non-linear 
fit procedure involving a free offset^^, of the form 



LK(Bo,T)^A 



Cm*T/Bo 
sinh (Cm* r /Bo) 



-Const. 



(A2) 



using three fit parameters. A, m* and Const. We will demon- 
strate in the following that such a procedure always leads to 
an overestimate of the mass m*. 

Effectively, adding a free offset to the fitting procedure is 
equivalent to the assumption that the power spectrum is equal 
to \d{F — Fo)\^ + \f{F)\^, i.e. the sum of the signal and noise 
in frequency space. However, as stated in eq. lAll it is not the 
case, but rather, we have that 



* I LK(T) + Noise | 

LK(T), m*/m = 12 

LK(T) fit, m*/m = 1 1 .4 + .6 

LK(T) + C fit, m*/m = 14 ± 2 . 
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FIG. 1 1 : Example of the two types of non-linear LK fits to the data. 
The absolute value of the sum of the LK function (blue curve) and 
random gaussian noise is shown in black stars, using a quasiparticle 
mass of 12 me. To this data was fitted the LK expression using a 
two parameter fit (red curve), which yielded a small underestimate 
of 11. 4± .6 ntg. A three parameter fit was also performed, featuring 
a free offset (green curve), which yielded an overestimate of the mass 
of 14±2m„. 



and the difference between the two sides of this inequality de- 
pends on the temperature. As the noise amplitude increases, 
the difference between these expressions increases, and the 
signal becomes 'buried' into the noise. Consequently, the 
width of the peak emerging from the noise at low temperature 
appears to decrease, yielding a larger mass m* when using a 
free offset. It is however not what happens when one uses the 
usual two parameter fit. 

In Fig. [TT] we illustrate this situation. We took the nor- 
mal LK distribution (blue curve), using a mass of 12 m^, to 
which we added random gaussian noise, and took the abso- 
lute value (the square root of the power spectrum), shown in 
black stars. We used 500 data points in order to show that 
the problem does not arise due to under-sampling. We then 
performed the two different types of fits, with (green curve) 
and without (red curve) a free offset. The resulting masses 
are different, where the two parameter fit obtains a slight un- 
derestimate (5%) while the three parameter fit overestimates 
the value by 17%. While neither of the two methods yields 
the correct value, the three parameter fit always deviates more 
strongly, and always overestimates the value. 

This effect depends strongly on the signal to noise ratio. 
Consequently, if the signal to noise ratio decreases gradually 
as a function of magnetic field, then the overestimates gener- 
ated by three parameter fits usually increases gradually, giving 
the impression of a field varying mass. In order to demonstrate 
the serious nature of this problem, we performed a simulation 
of typical dHvA data (similar to that of the a2 band), which 
featured beat patterns, a Dingle amplitude reduction factor and 
a constant level of uncorrected gaussian distributed noise of 
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FIG. 12: Simulations of dHvA data featuring a field dependent sig- 
nal to noise ratio, and its mass analysis using a three parameter fit 
featuring a free offset. The top panel presents noiseless simulated 
dHvA data, featuring two frequencies of 4.20 and 4.25 kT, leading 
to beat patterns, a mean free path of 1 /jm and a quasiparticle mass 
of 12 lUe. The following two panels show the same data with two 
different added noise levels. These sets of data were analysed with 
three parameter fits, the results of which are presented in the bottom 
three panels. 



LK{BJ)D{B) 



cos (271^1 /B )+ cos (27:^2 /B ) 



-/(B), (A4) 



where LK is the LK function, which uses a mass of 12 nie, D 
is the Dingle factor, which uses a mean free path of 1 /jm and 
a frequency of 4.2 kT and f{B) is the noise. 

In Fig. [12] top three panels, we show simulated raw data, 
the first without noise and the following two with different 
noise amplitudes. The bottom three panels present the respec- 
tive extracted field dependent masses when using three param- 
eter fits. The procedure for this calculation was the same as 
that which was used for the experimental data, described ear- 
lier. We found strong systematic errors in both cases where 
noise was added to the data, which are largest at each node of 
the beat pattern in the region where the overall signal to noise 
ratio is the lowest. We performed the same analysis while us- 
ing two parameter fits, not shown here, and this problem was 
not encountered. 



In Sr3Ru207 dHvA data, the signal to noise ratio decreases 
near the metamagnetic transition, mainly due to a combination 
of the Dingle reduction factor and scattering which increases 
near the QCEP and reduces the amplitude of quantum oscil- 
lations. Consequently, it is natural that this problem could 
have arisen specifically in this field range. However, it did not 
arise in the analysis of the data reported in this paper which 
featured a much higher signal to noise ratio, where the result- 
ing offset was generally close to zero. Such fits generated the 
same results as two parameter fits. 
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the critical point is approachedil At zero field and without impu- 
rity scattering, the temperature dependent resistivity contribution 
from the 72 pocket, which comprises a maximum of only one fifth 
of the Fermi surface, would be shorted out by the other pockets. 
However, at finite field as well as in the limit of residual scatter- 
ing dominating the resistivity a realistic calculation based on the 
multi-band Fermi surface of Sr3Ru207 shows that the transport 
could, in principle, be accounted for by a factor of 10 increase in 
the Y2 mass. 
Private communication. 



